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We have measured independently both s-wave and p- 
wave cross-dimensional thermalization rates for ultracold 40 K 
atoms held in a magnetic trap. These measurements re- 
veal that this fermionic isotope has a large positive s-wave 
triplet scattering length in addition to a low temperature p- 
wave shape resonance. We have observed directly the p-wave 
threshold law which, combined with the Fermi statistics, dra- 
matically suppresses elastic collision rates at low tempera- 
tures. In addition, we present initial evaporative cooling re- 
sults that make possible these collision measurements and are 
a precursor to achieving quantum degeneracy in this neutral, 
low-density Fermi system. 



While many examples of quantum degenerate 
fermionic systems are found in nature (electrons in met- 
als and nucleons in nuclear matter, for example), low 
density Fermi systems are exceedingly rare. However, 
techniques similar to those that led to the observation of 
Bose-Einstein condensation (BEC) in atomic systems 
can be exploited to realize a quantum degenerate, dilute 
gas of fermionic atoms. Novel phenomena predicted for 
this system include the suppression of inelastic collisions 
H , linewidth narrowing through suppression of sponta- 
neous emission ||, and the possibility of a phase transi- 
tion to a superfluid-like state at sufficiently low temper- 
atures Q . Just as was found in the case of BEC in alkali 
atoms, knowledge of the binary elastic collision cross- 
sections and interatomic potentials is crucial for experi- 
ments on fermionic species. Both evaporative cooling and 
the prospect of fermionic superfluidity depend on the cold 
collision parameters. For example, accurate prediction of 
magnetic- field Feshbach resonances ||, which could be 
used to realize Cooper pairing of fermionic atoms, hinges 
on detailed understanding of the interatomic potentials. 
In this letter, we present measurements of elastic collision 
cross sections for evaporatively cooled 40 K, including a 
direct observation of p-wave threshold behavior and the 
resultant strong suppression of the collision rate. 

Among the stable fermionic alkali atoms, 40 K yields 
the greatest range of possibilities for evaporative cool- 
ing strategies and interaction studies because of its large 
atomic spin, F (F=9/2 and F=7/2 hyperfine ground 
states). In addition to having a large number of spin 
states that can be held in the usual magnetic traps, 



potassium has two bosonic isotopes, 39 K and 41 K, which 
could be used for future studies of mixed boson-fermion 
dilute gases. Forced evaporative cooling which has 
proven essential for achieving quantum degeneracy in 
bosonic alkali gases, relies on elastic collisions for rether- 
malization of the trapped atomic gas. However, evapo- 
rative cooling strategies for fermionic samples are com- 
plicated by the fact that atomic collisions at these low 
temperatures (100's of /iK and below) are predominantly 
s-wave in character for bosonic atoms, while the Pauli ex- 
clusion principle prohibits s-wave collisions between spin- 
polarized fermions. Evaporative cooling for fermionic 
atoms must then proceed either through p-wave collisions 
or through sympathetic cooling ||. 

Since the s-wave and p-wave binary elastic cross sec- 
tions arc not well known for 40 K pf-pH, we have made 
measurements of elastic collision rates in a magnetic trap. 
The Fermi-Dirac quantum statistics of 40 K provide a 
unique opportunity to observe p-wave collisions directly. 
Exploiting this fact, we have seen evidence for a p-wave 
shape resonance and have observed threshold behavior 
of the p-wave cross section. To our knowledge, this mea- 
surement represents the first direct verification of the p- 
wave threshold law for neutral scatterers. 

We use a double-MOT (magneto-optical trap) appara- 
tus [|l2| to trap and pre-cool 40 K atoms prior to load- 
ing them into a purely magnetic trap. Operation of 
the MOT's employs two MOPA (master-oscillator power- 
amplifier) diode laser systems , each frequency stabi- 
lized using the DAVLL (dichroic absorption vapor laser 
lock) technique [|l4| which accomplishes a large frequency 
range for locking. Also, we have developed an atom 
source [H| that is enriched with 5% 40 K (whose natu- 
ral abundance is 0.01%). This system allows us to trap 
any of the three potassium isotopes, and to trap 10 s 40 K 
atoms (four orders of magnitude more 40 K atoms than 
previous efforts [[16)). 

Immediately before loading the sample into the mag- 
netic trap, the atoms are cooled to approximately 150 fiK 
during a Dopplcr cooling stage of the MOT in which the 
trapping light is jumped closer to resonance. Further, 
an optical pumping pulse transfers the majority of the 
atoms into magnetically trappable states in the F=9/2 
hyperfine ground state. The atoms are then loaded into 
a cloverleaf magnetic trap [ p"7| and after an initial evap- 
orative cooling stage, we are left with roughly 10 7 40 K 
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atoms at 60 /iK and a peak density of 10 9 cm~ 3 . The 
cloverleaf magnetic trap provides a cylindrically symmet- 
ric harmonic potential, with a characteristic radial fre- 
quency v r — 44 ± 1 and an axial frequency v z of 19 ± 1 
Hz for loading. The lifetime of the atoms in the mag- 
netic trap, limited by collisions with room-temperature 
background atoms, has an exponential time constant of 
300 ±50 s, giving ample time for thermal relaxation stud- 
ies as well as for evaporation. The radial frequency and 
the bias magnetic field can be altered smoothly by chang- 
ing the current through a pair of Helmholtz bias coils. 

We determine elastic collision cross sections from mea- 
surements of cross-dimensional thermalization rates fl§| l 
in the magnetic trap. The sample is taken out of thermal 
equilibrium by changing v r through a ramp of the bias 
coil current. For the measurements reported here v r lies 
between 44 and 133 Hz. The change in v T occurs adiabat- 
ically (slow compared to the atomic motion in the trap) 
but much faster then the rate of collisions between atoms. 
Since the axial frequency is essentially unchanged, energy 
is added to (or removed from) the cloud in only the radial 
dimension. Elastic collisions then move energy between 
the radial and the axial dimensions, and the thermal re- 
laxation is observed by monitoring the time evolution of 
the cloud's aspect ratio. 

To avoid perturbations to the image due to the spa- 
tially dependent magnetic fields, the trap is turned off 
suddenly and the cloud is imaged after 2.7 ms of free ex- 
pansion. The aspect ratio of the cloud is observed via 
absorption imaging using a 9.1 fis pulse of light resonant 
with the 4Si /2 , F=9/2 to 4P 3/2 , F=ll/2 transition. Op- 
tical depth is calculated from the image captured on a 
CCD array and then surface fit to a gaussian distribution 
to find the rms cloud size in both the radial and axial di- 
mensions. An example of the cloud evolution following 
a change in trap potential is shown in Fig. 1. Since the 
expanded cloud sizes are proportional to the square root 
of the cloud energy in each dimension, the exponential 
time constant for the redistribution of energy, r, can be 
extracted from an appropriate fit to the aspect ratio vs 
time. To rule out significant relaxation through trap an- 
harmonicities, we have verified that the relaxation rate 
1/r scales linearly with the number of trapped atoms N. 

To obtain the elastic collision cross section a from our 
measurements of thermal relaxation rates, we use the re- 
lation: 1/r = -^ncrv, where n is the density- weighted 
density of the trapped atoms given by ■h J n(r) 2 d 3 r, v is 
the rms relative velocity between two atoms in the trap, 
and a is the calculated average number of binary colli- 
sions per atom required for thermalization. The product 
nv depends on both the size and temperature, T, of the 
trapped sample. These are measured by observing the 
expansion of an equilibrated sample after release from 
the magnetic trap. The rate of expansion yields the tem- 
perature, while an extrapolation back to the release time 
gives the initial sizes. Using the trap potential calculated 



from the field coil geometry we have checked that the 
measured initial sizes and temperatures are consistent to 
within their uncertainties. 

The mean number of collisions each atom undergoes, 
a, during one relaxation time constant was determined 
from a numerical simulation of the experiment using clas- 
sical Monte Carlo methods [[l9| . For a harmonic trapping 
potential, the relaxation simulation yields a s — 2.5 for 
s-wave collisions and a p = 4.1 for p-wave collisions. The 
ratio a p /a s can also be determined analytically through 
an integration over the angular dependence of scattering. 
This gives a p /a s = 5/3, consistent with the Monte Carlo 
results. 

The primary results of this paper are shown in Fig. 2. 
While ordinarily one cannot measure higher order partial 
wave contributions to the collision cross section directly, 
the Fermi-Dirac statistics of 40 K allow us to probe p- 
wave and s-wave interactions independently. The p-wave 
cross section a p is determined from measurements using 
a spin-polarized sample (\F = 9/2, mp = 9/2 > atoms), 
where s-wave collisions are prohibited by the quantum 
statistics. The s-wave cross sections <r s are determined 
from data obtained using a mixture of two spin states, 
|9/2, 9/2 > and |9/2, 7/2 >. The magnitude of the p- 
wave cross section is surprisingly large and we find that 
40 K has a p-wave shape resonance at a collision energy of 
roughly 280 [iK. At temperatures well below the resonant 
energy (less than 30 ^K), a fit to <r p vs T gives a p cx 
T 2 0±0 ' 3 . Thus, we have directly observed the expected 
threshold behavior a p oc E 2 . In contrast, a s exhibits 
little temperature dependence. With these very different 
temperature dependencies, the collision rate changes by 
over two orders of magnitude at our lowest temperatures 
depending on the spin mixture of the fermionic atom gas. 

To explore this effect further we measure the ther- 
malization rate vs spin polarization at 9 fiK (see Fig. 
3). We control the relative populations of |9/2, 9/2 > 
and 1 9/ 2, 7/2 > atoms in a two-component cloud with 
a microwave field that drives transitions to untrapped 
spin states in the F=7/2 ground state manifold. The 
trap bias magnetic field breaks the degeneracy of the hy- 
pcrfinc ground-state splitting (1.286 GHz at zero field 
ppf ) so that the different spin-states can be removed 
selectively (see Fig. 3 inset). For the data shown in 
Fig. 3 the fraction of atoms in the |9/2, 9/2 > state 

f mF =9/2 = N mF=9 "T+N2 F= 7/2 ' was varied smoothly from 
70 to 100% by varying the power of an applied microwave 
field (frequency swept) that removes a portion of the 
9/2, 7/2 > atoms. 

The thermalization of mixed spin-state samples de- 
pends on both s-wave and p-wave collisions. The data 
in Fig. 3 can be fit to a simple model given by: 

2 2 C 2 
1/r = Ui, 2 ( — cr s H ^-)v + (m,i + n 2 , 2 ) — v P v, 

GL q C^p " ^*p 
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where riij is the density-weighted density between two 
species given by n^j — N \ N J ni(r)rij(r)d r and the 
subscripts 1 and 2 stand for the two relevant spin states. 
Since the magnetic moments of the |9/2, 9/2 > and 
1 9/2 , 7/2 > atoms are only slightly different we make the 
simplifying assumption that these states have identical 
spatial profiles in the trap. A fit using the above model 
with a s and the ratio a p /a s as free parameters shows 
good agreement with the data in Fig. 3. In addition 
to demonstrating the type of control over collision rates 
that is available in a trapped gas of fermionic atoms, this 
measurement of a p /a s at low temperature provides a sen- 
sitive constraint on the triplet scattering length. The s- 
wave cross-sections shown in Fig. 2 were extracted using 
the above equation, however at these low temperatures 
Op is relatively small and the measured thermalization 
rates are due primarily to s-wave interactions. 

To compute the scattering cross sections for com- 
parison with these data, we first identify singlet and 
triplet potassium potentials that have been deter- 
mined spectroscopically. At large interatomic separation 
R these potentials are matched smoothly to the long- 
range form of Marinescu et al. p2l . We add an additional 
correction to these potentials' inner walls fTl]] , enabling 
us to tune the scattering lengths over their entire ranges 
-co < a < oo. We set the singlet scattering length's 
value at a s — 104ao Jll[] where a is the Bohr radius, but 
leave the triplet scattering length at as a free parameter 
to be determined by the experiment. We note that the 
present experiment is relatively insensitive to the value 
of a s since even the |9/2, 9/2 > +|9/2, 7/2 > process is 
strongly triplet-dominated and no singlet resonance oc- 
curs near threshold; indeed, varying a s over its range of 
uncertainty, 101 < a s < 107 |0], does not change the fit. 

After computing cross sections as a function of colli- 
sion energy, we determine temperature-dependent cross 
sections by computing a thermal average over collision 
events, weighted by the collision energy. Using this 
type of thermal averaging to account for a temperature- 
dependent cross section is supported by Monte Carlo 
studies p9| . To make a fit to the data, we compute x 2 
while floating both at and a multiplicative factor e which 
scales simultaneously the computed cr s (T) and a p {T). 
This factor is required to accommodate a ±50% system- 
atic uncertainty in the experimental determination of ab- 
solute cross sections (primarily from N). Our global best 
fit occurs for at = 157±20ao and e — 1.6, with a reduced 
X 2 of 3.8; the corresponding cross sections are plotted as 
lines in Fig. 2. The uncertainty in at reflects a doubling 
of the fit x 2 an d includes a ~ 2ao uncertainty arising 
from varying Cq over its range 3600 < Cq < 4000 a.u. 
p2| . Our nominal potential gives a p-wave shape reso- 
nance at ~ 280/iK in collision energy, with an asymmetric 
lineshape whose FWHM is ~ 400^iK. 

The relatively small uncertainty on the value of at is 
attributable to the fact that we can simultaneously fit 



s-wave and p-wave collision data having little relative 
uncertainty. The value of at for 40 K determined here 
does not agree well with reference 113], highlighting the 
importance of low temperature data in determining ac- 
curate potentials. Our measurement is however in good 
agreement with the value at = 194+42 2 obtained recently 
from an analysis of photoassociation spectroscopy of 39 K 
|ll| |. This agreement between two fundamentally differ- 
ent experiments is very encouraging, and suggests that 
the potassium scattering lengths are now fairly well de- 
termined p3[ . Using our new value, we have tabulated 
in Table 1 the resulting triplet scattering lengths for col- 
lisions between the different potassium isotopes. 

One of the important applications of these thermaliza- 
tion rate measurements is in determining the feasibility of 
various evaporative cooling strategies for this fermionic 
atom gas. The large p-wave cross section makes evap- 
oration of a spin-polarized sample possible, but only at 
T > 20fiK. To reach lower temperatures, which is likely 
to be necessary for producing quantum degenerate sam- 
ples, evaporation of a mixed spin-state sample should 
work well given the suitably large s-wave elastic collision 
cross section. 

Indeed, for the data presented in this letter we vary the 
temperature of the sample through forced evaporation. 
To facilitate evaporative cooling, we increase the colli- 
sion rate by ramping to a v r = 133 ± 5 Hz trap. Evap- 
oration then proceeds using a mixed spin state sample 
and applying a microwave field to remove selectively the 
most energetic atoms (in both spin states). With these 
initial attempts at evaporatively cooling of K we can 
lower the sample temperature from 100 [iK to 5 fiK. We 
see runaway evaporation, where the collision rate in the 
gas increases as the temperature decreases. While the 
samples discussed in this work are still far from quantum 
degeneracy, this initial success and our measurement of 
a large s-wave elastic collision cross section bode well for 
future progress toward this goal. 
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FIG. 1. Example of thermalization data. The inset shows 
the axial size, z rms (•), and the radial size, x rms (o), imaged 
after 2.7 ms of free expansion, relaxing as the trapped atoms 
rethermalize via elastic collisions. At time=0 the cloud is 
taken out of equilibrium by changing u r from 133 to 44Hz. 
A fit (line) to the aspect ratio vs time is used to extract the 
relaxation rate. 



4 




FIG. 2. Elastic cross sections vs. temperature. The s-wave 
cross section (o), measured using a mixture of spin states, 
shows little temperature dependence. However, the p-wave 
cross section (•), measured using spin-polarized atoms, ex- 
hibits the expected threshold behavior and is seen to vary by 
over two orders of magnitude. The lines are a fit to the data, 
as described in the text, yielding at = 157 ± 20ao. 



FIG. 3. Dependence of a on spin composition at T — 9fiK. 
A quantity proportional to a, 1/tN, is measured as a func- 
tion of the fraction of atoms in one of two trapped Zeeman 
spin states, fm F =9/2- The inset shows the number of atoms 
remaining after application of microwaves at the indicated fre- 
quency. The three features correspond to removal of trapped 
atoms in particular spin states and can be used to measure 
or control the spin composition of the atom gas. 



TABLE I. Triplet scattering lengths a t in Bohr radii for 



collisions between potassium isotopes. 



Isotopes 


Nominal at 


Range 


40 + 40 


157 


136 <a t < 176 


39 + 39 


-44 


-80 <a t < -28 


41+41 


57 


49 < a t < 62 


39 + 40 

39 + 41 

40 + 41 


3600 
164 
93 


a t > 500 or a t < -900 
140 <a t < 185 
83 < a t < 99 
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